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At second order in gradients, conformal relativistic hydrodynamics depends on the viscosity r\ 
and on five additional "second-order" hydrodynamical coefficients rn, K, Ai, A2, and A3. We derive 
Kubo relations for these coefficients, relating them to equilibrium, fully retarded 3-point correlation 
functions of the stress tensor. We show that the coefficient A3 can be evaluated directly by Euclidean 
means and does not in general vanish. 
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Results from the RHIC experiments, particularly the 
measurement of a large transverse flow 1], appear to 
show that the Quark-Gluon plasma can be well described 
by hydrodynamics with a surprisingly small viscosity [2| . 
A major future goal for heavy ion experiment and the- 
ory is to quantify how small the viscosity of the plasma 
is. This requires the numerical treatment of relativistic 
viscous hydrodynamics [3J . It has long been known [4]-|6[ 
that the relativistic Navier-Stokes equations are acausal 
and unstable. But the Navier-Stokes equations are just 
the result of a first-order (Chapman- Enskog [7]) expan- 
sion in gradients. Extending the expansion to second or- 
der yields numerically stable equations after a certain re- 
organization is applied |5| . The drawback is that it adds 
unknown coefficients. In the conformal case (which we 
will consider for simplicity) , besides the equation of state 
P(e) at zero order and the shear viscosity n at first order, 
there are five new transport coefficients: rn, K, Ai, A2, A3 
in the notation of [8]. These have been evaluated in 
strongly coupled A/"=4 super- Yang-Mills theory in the 
limit of many colors [8j, [9| and at leading order in weakly 
coupled QCD [lfjj. In each case A3 = at lowest order 
in the respective (strong or weak coupling) expansion. 

Baicr et al have also presented Kubo formulae for 
two of these coefficients, tu and k, which relate them 
to well denned, equilibrium correlation functions of the 
stress tensor. Presumably, the remaining three coeffi- 
cients Ai^.3 can also be expressed in terms of stress tensor 
correlation functions. Doing so would put the definition 
of these coefficients on a solid footing and might aid in 
their physical interpretation and their theoretical calcu- 
lation. In the remainder of this paper we will derive such 
Kubo relations for the three remaining second-order co- 
efficients. We do this first by showing how the first and 
second order hydrodynamic coefficients can be related to 
the stress tensor in a background spacetime with pertur- 
batively small geometrical curvature. Then we expand 
in the metric as an external background field a la Kubo 
11] and derive a relation between Ai^.3 and certain fully 
retarded 3-point stress-tensor correlation functions. This 
allows us to determine the previously unknown pertur- 
bative behavior of the coefficient A3 (which is not zero) 
and to say something about its physical interpretation. 

We restrict attention to conformal fluids mostly to sim- 



plify the presentation; in the nonconformal case there 
are more coefficients [12 ] but there are no conceptual or 
technical obstacles to treating this case with the same 
methodology developed here. 



CONSTITUTIVE RELATIONS FOR SECOND 
ORDER COEFFICIENTS 



We begin by defining the second order coefficients. The 
expectation value of the stress-energy tensor operator for 
a fluid can be decomposed in terms of a local equilibrium 
piece and an extra piece, 
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Here g M „, e, P, iP are the spacetime metric (in the mostly- 
plus convention), energy density, pressure as given by the 
equation of state, and flow 4-velocity. We work in the 
Landau-Lifshitz frame, u M IP" = 0, which makes the di- 
vision between T^i" and IP" unique; we normalize tp so 
that u^tP — — 1. While iP, e, and g^ u are ordinary func- 
tions of x, T M " is a Heisenberg-picture operator; (T M ") 
represents its trace in the density matrix describing the 
fluid. 

The key idea of hydrodynamics is that, for a system 
which varies slowly in space and time, IP" arises only due 
to the nonuniformity of the system and should therefore 
be expressible in terms of a gradient expansion in that 
nonuniformity. To write out IP" to second order, we in- 
troduce some notation. We define A M " = <p" + iPu" ', 
which is the projector to spatial directions in the local 
rest frame. Angular brackets around a pair of Lorentz in- 
dices, ^"^ , mean that the indices are to be symmetrized, 
space-projected, and trace-subtracted; that is, 
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The shear and vorticity tensors are defined as 

a p,v s 2\J^u v \ (3) 

Or = lA^A^V^-V^). (4) 

W" and B} LVa $ are the Ricci tensor and curvature tensor 
respectively. In terms of these quantities, the most gen- 
eral form for IP" compatible with conformal symmetry 
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EXPANSION IN BACKGROUND GEOMETRY 

We derive Kubo relations for Ai etc. by considering 
a system where some nonuniformity, either in the ini- 
tial conditions or in the spacetime geometry, forces a 11 " 
etc. to be nonzero. It is particularly convenient to con- 
sider an initially uniform, equilibrium system in flat space 
but to introduce perturbatively weak and slowly vary- 
ing spacetime nonuniformity which causes the fluid to 
experience shear and vorticity. Writing the metric as 
9nu{x) = rj^ v + h lliy (x) (rj^ the flat-space metric), one 
expands perturbatively in /i M „. Since h^ v couples to the 
stress tensor T^ v , this generates an expansion in cor- 
relation functions of multiple stress tensors, whose co- 
efficients are the response of the stress tensor to fluid 
nonuniformitics. 

Consider the expectation value (T A " y (0)) for a system 
initially (time to <C 0) in equilibrium at temperature T, 
subject to a spacetime dependent metric perturbation 
h a p{x), with hfj,i,(t < to) = 0. The stress tensor is deter- 
mined by 

(T^(0)) = Tr e~ fiB f cxp ( f dt'iH[h(t')]) T" v 



xTexpl / dt" {-i)H[h{t")\\ (6) 

(with Texp and Texp the anti-time ordered and time- 
ordered exponentials respectively, H[h(t)] the Hamil- 
tonian, showing explicitly its dependence on the met- 
ric, and (3—T^ 1 the inverse temperature). This is best 
treated using the Schwinger-Keldysh (closed time path) 
formalism (see 13|, [l4[; we follow the conventions in 
|14j). We introduce independent metric perturbations 
for the T-ordered and T-ordered evolution operators in 
the above expression and define the generating functional 



W[hxM] = lnTre-^fexpN f dt'H[h 2 (t')] 
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One then defines the average metric perturbation h r = 
hl + h2 and stress tensor T r = Tl ~^ T2 , and the difference 



variables h a = hi — hi, T a = T\ 
respect to h a gives T ri explicitly 
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We use such a variation to pull down the T M "(0) factor 
we want to evaluate. After taking this h a derivative, we 
set h a = and h r = h, since we are interested in the case 
of a classical background value hi = h% = h r = h. (The 
difference h a represents possible quantum fluctuations in 
the metric which we do not want to consider.) We then 
expand order by order in h rfiu to obtain a series expansion 
of {T^ u ) in powers of h. Explicitly, we find 

(TDh = GT(0)-i/ ' d*xG%' a f l (p,x)h a f,(x) (9) 

+ i J d^xd% GZa aP ' lS (0, x, y)h a p(x)h lS (y) 

plus terms of order h 3 . Here G£%'""'"'(p,x,...) is the 
correlation function of one T r and or more T a 's, 
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The expectation value is with respect to the flat-space, 
equilibrium density matrix. G ra ,,, is a fully retarded 
correlation function [14| . which is a nested commutator, 
from earliest to latest time, with T r at the last time and 
innermost in the commutator, eg when x° < y° < . . . < 
the correlator is ([T{x), [T(y), [. . .T(0)]]]). Here (c.t.) 
refers to the contact terms which are built into our def- 
inition of the n-point stress tensor correlation functions. 
This is discussed in [15|; the contact terms turn out not 
to be important for evaluating r\ but they will contribute 
to the evaluation of Ai o,3- 



KUBO FORMULAE 

First we review the derivation of Kubo formulae for 
the "linear" transport coefficients rj,Tn_,K [8J. Consider 
(T xy ) in the presence of h xy (z,t). According to Eq. ©, 
at first order 
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Using Eq. ((T|) and V M T M,y = (energy-momentum con- 
servation), we derive that u l = at 0(h). We then 
evaluate a xy , u ■ Ver'^ etc. explicitly for this u M and 
h^v, finding for instance that a xy = dth xy . Substituting 
into Eq. ([5]), we find 
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defining G xy - xy (oj,k) = J dW^'-^G^'^CO, -x) and 
equating Eqs. (J11I12I) order by order in derivatives, we 
find 
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(d 2 u G x r y ' xy (w, k) - dl G xy ' xy ( U , k)) l =Q=k .(15) 



These reproduce the Kubo relations obtained by [8[. 

To obtain higher order Kubo formulae for the nonlinear 
coefficients, we continue this procedure to 0(h 2 ), for a 
background choice which allows nonzero shear flow and 
vorticity. To do so, we will consider Tl xy arising when 
h X z(i), h yz (t), h x o(z), and h y o(z) are nonvanishing[19]. 
By our choice of indices, a xz = d t h xz , a yz = d t h yz , fl xz = 
—d z h x0 /2 and Vt yz — —d z h y o/2 arise at 0(h), but H xy 
will automatically only arise at order in h 2 , and will not 
receive contributions at this order from the (e+P)u x u y 
term. Since g xy is explicitly 0(h 2 ), we also only need the 
equilibrium value of P[20(. Explicitly evaluating Eq. ([5]) 
in this background to second order, we find 

(T xy ) = P(h xz h yz - h x0 h y0 ) +r](h xz h yZit + h XZit h yz ) 
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Equating with the h 2 part of Eq. ([9]) , and defining 

(17) 
we find the following Kubo relations: 
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These Kubo relations are our main result. 

We also find extra Kubo relations for 77, k, and t w : 
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G x r y a 'a°' v0 (P,q), (22) 
G^'.v-Cp.ff), (23) 



where besides the differentiated variable all other p, q 
components are taken to zero first. These extra relations 



require inter-relations between G^ ,a/3 and G^"^' 7 " 5 . 
Each extra Kubo formula involves one stress tensor at 
zero external 4-momentum, arising from an undifferenti- 
ated h^ in Eq. (fll)|). We can always force h^ v = at 
x = where T xy is evaluated by a coordinate "gauge" 
choice. The invariance of the theory to such gauge 
choice enforces (Ward) relations between two point func- 
tions and three point functions with a T£ v at zero 4- 
momentum. Consider a stress tensor two-point function 
in a spacetime- independent, background h^: 



(Tr(0)T^(x)) h = iG%>'#(p,x) (24) 
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The gauge change which eliminates h^ v is x M — > a;^ + £ M 
with ^ )W + £„ !jU = hfj, v . Applying the gauge change to 
the lefthand side of Eq. (|24|) , we re-express it in terms of 
h^v and the flat-space correlation functions; 

M(['T , ^(p) + 0'«'')] 

+ [rT &%**>$) + (a o /?)] +(70 *)' 
= 2/ M G^ 3 ^(p,0). (25) 

Choosing /if = xy, a/3 — xz, j8 = yz, 

G x r y a xy (p) + G X T 7 \V) = 2G xy aa x ^ yz (P, 0) • (26) 

Now d u G xy > xy = irj by Eq. JT3) and d u G xz ' xz = ir] by 
rotational invariance, soEq. (l2lT) follows. The same pro- 
cedure applies for the other linear coefficients. 

DISCUSSION 

Our derivation had two goals. First, we wanted re- 
lations, shown in Eqs. (|18I19I20[) . for the second-order 
nonlinear transport coefficients in terms of equilibrium 
energy-momentum tensors. Second, we hoped that these 
relations would shed some light on the nature or proper- 
ties of these transport coefficients. The most mysterious 
of these transport coefficients is A3, which is found to 
vanish in 7V=4 SYM theory in the limit of many colors 
and large coupling [9j and which is zero at order g" 8 in 
the weak coupling expansion, the order where Ai^ are 



nonzero [10|. Is it identically zero? Romatschke [12 1 
studied this problem (among others) using a generalized 
entropy current and showed that A3 is related to a cer- 
tain modification of the entropy density in the presence 
of vorticity. Our Kubo relation allows for a direct evalu- 
ation of A3 in weakly coupled field theory. 

k and A3 have expressions involving space but not 
time derivatives of stress tensor correlation functions. 
We may immediately set u = in Eq. (fblj) and 
p°,q° — in Eq. ([2"0")l . The frequency-domain fully 
retarded function G ra ... is the analytic continuation of 



Ge the Euclidean correlation function. In particular, 
^^■^(-iwi.-wa) = i^G^'^iu!,^) for Mat- 
subara frequencies Wi,2 = 27rTni ) 2- Here no is the num- 
ber of indices /j,,v,a,/3,cr,T which are 0, since there is 
a factor of i arising from the Euclidean continuation of 
a index. This relation shows that the zero-frequency 
raa and Euclidean correlation functions are equal up to 
factors of i . Hence 



A3 = 4 lim 



f>,q^Q dp z dq z 
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One usually considers such Euclidean correlation func- 
tions to carry only thermodynamical information; A3 
should not be thought of as a dynamical coefficient but 
as a thermodynamic response to vorticity[21]. 

At weak coupling we can directly evaluate Eq. (j2"7) di- 
agrammatically in the Matsubara formalism. This con- 
trasts with the case of n,rn,Ai, and A2, where time 
derivatives mean that G raa must be evaluated at small 
nonzero frequency where the continuation cannot be so 
simply applied. Therefore the weak coupling expansion 
of A3 (and k |15j |) will start at g°, while the expansions 



17] 



for A12 can involve inverse powers of g [1 

We have evaluated the correlation function in Eq. (j2"7| 
for a one-component scalar field theory at leading or- 
der in weak coupling. Two diagrams contribute; a tri- 
angle diagram, d Pz d qm (T^(-p - q)T x0 (p)^ (q)) = ^ 
and a contact term involving X x0v ° = 2dT x0 jdg. 
d Pz 8 qz (T x y(-p - q)X x °y°(p + «)) = §. Hence 



8/0) 



A3 
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1 weak-coupled real scalar field. 



(28) 



We get the same answer using the scalar field stress tensor 
from 18|. The important observations are that A3 can 



be quite easily evaluated at weak coupling via Euclidean 
techniques, and the result is not in general zero. 
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